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We have developed an experimental method to modify the single-particle dispersion using peri¬ 
odic modulation of Raman beams which couple two spin-states of an ultracold atomic gas. The 
modulation introduces a new coupling between Raman-induced spin-orbit-coupled dressed bands, 
creating a second generation of dressed-state eigenlevels that feature both a novel Stt spin-orbit 
coupling and a pair of avoided crossings, which is used to realize an atomic interferometer. The 
spin polarization and energies of these eigenlevels are characterized by studying the transport of a 
Bose-Einstein condensate in this system, including observing a Stueckelberg interference. 
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In ultracold atoms, laser-induced synthetic gauge fields 
have been used to modify the single-particle dispersion 
and realize physics such as synthetic magnetic fields [T], 
spin-orbit (SO) coupling [5] and the Spin-Hall effect [3], 
the Hofstadter/Harper and Haldane Hamiltonians [IHS]. 
Raman-coupling between spin states of ultracold atoms 
which modifies the free-particle dispersion relation [7] has 
been the basis of many of these works Eiin], and has re¬ 
sulted in a rich field of studies in Rashba-Dresselhaus 
type SO coupling for both Bose-Einstein condensates 
(BECs) and degenerate Fermi gases pTiHTH] . 

In this work, we investigate adding a time-dependent 
modulation to the intensity of the Raman-coupling as a 
new way to engineer the dispersion relations of ultracold 
atoms. We demonstrate two results to show the effec¬ 
tiveness of this technique: (1) by observing a Stt rotation 
of the spin polarization in the spin-momentum locked 
dressed eigenlevels and dispersion relation, which we call 
a Stt SO coupling, and (2) using the eigenlevels in the 
dressed state picture of the modulated Raman system to 
create an atom interferometer. 

Our experiment starts with using SO coupling of equal 
parts Rashba and Dresselhaus that is created by Ra¬ 
man coupling two spin states of an atom [2], shown in 
Fig. 0 (a) with a similar experimental setup to our pre¬ 
vious work (see Ref. |I7j). This SO coupling consists 
of two eigenlevels, Eu{q) and E^^q) for the upper and 
lower, both of which possess a spin-momentum “lock¬ 
ing” in quasimomentum, hq, space. We then create an 
inter-eigenlevel coupling by modulating the strength of 


the Raman coupling: UlR(t) = Dq -I- UImcos(2'k fmodt), 
where fmod is the modulation frequency, Oq is the un¬ 
modulated Raman coupling, and UIm is the modulation 
amplitude. 

In the theoretical analysis of this system, we first calcu¬ 
late the Rashba-Dresselhaus SO coupled eigenlevels from 
the Hamiltonian: 




Do/2 


Do/2 1^(5 _fc,)2 +5/2 


( 1 ) 


to find the SO eigenlevel structure, as pictured in Fig. 
(a) where 5 is the Raman laser detuning, m is the atomic 
mass, kr is the recoil momentum, and h is the reduced 
Planck’s constant {h = 2^1%). The third spin state is 
neglected due to the quadratic Zeeman shift. We define 
total atom spin polarization as 5 = {N^ — N^)/{N^+N^), 
where is the number of spin up (down) atoms in 

the BEG. 

Using the dressed state picture, these original SO 
eigenlevel structure {Eu{q) and EL{q)) forms the basis of 
the manifold which is replicated in every hf^od to form 
a set Ejj/R^n = Ejj/R+nhfmod where n = 0, ±I,±2,... in 
the modulated system. Considering just nearest neigh¬ 
bor eigenlevels (n = ±1) with a modulation frequency 
slightly larger than the eigenlevel energy difference near 
q = 0 {hfmod Eu{q Ks 0) - EL{q 0)), the ground 
eigenlevel EL^n=o has a pair of crossings with Eu^n=-i 
(the precise value of fmod determines the locations of 
the crossings in quasi-momentum space). We use these 
eigenlevels as the new bare states coupled by the modula¬ 
tion field with an effective coupling strength denoted Dq-. 
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FIG. 1. (a) The standard Rashba-Dresselhaus SO eigenlevels created via Raman coupling of two atomic mp ground states, here 
with ®^Rb and flo = l.SEr, S — 0. The dashed line shows the upper eigenlevel Eu shifted down by the modulation frequency 
fmod = 10.56 kHz. (b) Dressed-state diagram of the Stt SO-coupled eigenlevels created via modulation of the Raman coupling 
strength, which feature two avoided crossings (labeled A and B) here with the same parameters as (a) and Q.c = 0.27Er. The 
modulation-dressed upper eigenlevel is denoted E+, and the lower E-, with Rashba-Dresselhaus eigenlevels shown by dashed 
yellow lines. The green bars indicated the avoided crossings used as beam splitters (in this case approximately 50/50) in the 
interference experiments. In (a) and (b) the blue and red colors represent spin down {mp = 0) and up {mp = —1) respectively, 
(c) Measurements of the spin-polarization of a BEG for transport in both the un-modulated system shown in (a) and the 
modulated system in (b). The BEG starts at qi, labeled in (b), and falls under gravity {ap — 1680 kr/s) along the —q direction 
and nearly adiabatically follows the lowest energy eigenlevel. With no modulation, the full tt spin rotation is observed (black 
crosses). With modulation of flpi = l.SEr, a nearly Stt rotation is observed (purple circles) as the quasimomentum, hq, of the 
BEG goes from +hkr to —hkr. Solid black (purple dashed) lines are the calculated spin polarization of Ep {E-) given the same 
parameters as in experiment and Dc = 0.58i?r (calculated eigenlevel structure for this experiment is shown in Fig.|^. For this 
and following figures, a representative error bar indicates an average of 10% uncertainty in atom population in each spin due 
to technical noise. 


Solving at each q, we obtain the new modulation-dressed 
eigenlevels, labeled E+{q) and E-{q), as eigenlevels of 
the Hamiltonian 


( EL,n=q{.q) f^c/2 \ 

\ Hc/2 Eu,n=-l{q) ) 


( 2 ) 


These eigenlevels feature a double avoided crossing with 
gap size He. Both E+{q) and E-{q) feature a Stt rotation 
of the atom-state spin polarization as the quasimomen¬ 
tum goes from -\-hq to —hq, in contrast to the Itt rotation 
in the original Ejj and Ep. Transport of a BEC in the 
new eigenlevels E^{q) and E_{q) are basis of our exper¬ 
imental investigations and are used to observe both Stt 
rotation of the spin and to engineer an atom interferom¬ 
eter. 

For our experiments, we utilize transport methods 
developed in our earlier work m to study these 
modulation-dressed spin-orbit eigenlevels. Briefly, the 
initial state of the ®’^Rb BEC is prepared by starting with 
a spin-pure BEC in the mp = 0 state, and adiabatically 
turning on the Raman coupling to a fixed value Hq. This 
loads the atoms at the initial value qi « Ikp. The modu¬ 
lation of the Raman beams is then turned on, the optical 


trap holding the BEC is turned off, and then a force is ap¬ 
plied in the —q direction to accelerate the BEC at an av¬ 
erage rate ap through both the avoided crossings, labeled 
‘A’ and ‘B’ in Fig[^ The probability of transition to the 
upper Ep eigenlevel is determined by the Landau-Zener 
(LZ) probability Ppz = exp [—27r(Hc/2)^/(/ia/3)] where 
a = \dq/dt\ is the rate of acceleration at the avoided 
crossing and (3 is the difference of the slopes of the bare 
state energy levels m- The BEC is imaged after 15 ms 
time-of-flight expansion with a Stern-Gerlach field ap¬ 
plied to measure both the spin and momentum of the 
BEC. The ||) (|t)) state corresponds to ®^Rb hyperfine 
state \F = l,mp = 0(—1)). The recoil energy from the 
Raman lasers is Epi = h'^k^/2m = h x 3.68 kHz. 

Fig. 0 (c) shows the experimentally measured spin po¬ 
larization of the BEC accelerated by the gravitational 
force {ap = 1680 kr/s) along the ground dressed eigen¬ 
level. For Hm = 0, the measured BEC spin polarization 
(black crosses) follows the calculated spin polarization of 
Epiq) (black line). For the case of strong modulation 
(Hm = 1.3Er) the BEC nearly adiabatically follows the 
calculated spin polarization of E-{q). The Itt spin rota¬ 
tion of Ep is evident in the data with no modulation, and 
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a Stt spin rotation is observed in the data with modula¬ 
tion. The measured spin polarization of the BEC doesn’t 
perfectly match the calculated E_ eigenlevel spin polar¬ 
ization, as seen in the data after about 1 ms, due to im¬ 
perfect loading into the dressed state and non-adiabatic 
inter-eigenlevel transitions m- Nonetheless, this exper¬ 
iment demonstrates the viability of modulated-Raman 
coupling as a means to create a new type of SO coupling 
with a Stt spin rotation, different from the previously 
studied, static Raman-induced SO with Itt spin rotation. 

In addition to realizing a new type of Stt SO coupling, 
we use these modulated light induced synthetic gauge 
fields to engineer an atom-interferometer. Stueckelberg 
interference results when a quantum state is coherently 
split by LZ transition at an avoided-crossing, each split 
state travels along a different path, and then the states 
are recombined after another LZ transition at the same 
or a different avoided crossing [niin]. The final state 
depends on the energy difference of the two paths and the 
time it takes to traverse the path (together determining 
the total phase difference), and fringes in measurements 
of the final state will occur if the path or time is varied. 
The eigenlevels created using a traditional Raman cou¬ 
pling of internal atomic spin states, e.g. Fig.[^a), do not 
give a usable pair of avoided crossings [21]. However, in 
the modulated system investigated here, such a pair is 
readily realized, and an example is indicated by vertical 
lines A and B in Fig. [^b). 

For these experiments, the energy difference is that 
between the upper and lower eigenlevels and 

i?_, and the BEC splits first at quasimomentum 
QA and then is recombined at qb with both pro¬ 
cesses involving Landua-Zener transitions {Plz{oi) = 
exp [—27r(Hc/2)^/(/ia/3)]). The rate of acceleration, a, 
determines the time it takes to get from qA to gs, and 
the total phase difference acquired is calculated <i>(a) = 
[E+{q) — E-{q)] dq/{ha). For the above and follow¬ 
ing equations, the dependence on Hq, ftc, <5, and fmod 
has been suppressed for notational clarity. The output 
spin polarization, 5, for an atom moving through this 
interferometer with acceleration a is calculated as (see 
Appendix A): 

5(a) = 4 {Plz{o) - PLz{af) cos (^(a))-(l-2PLz(o:)f 

(3) 

We measured Stueckelberg interference fringes using 
eigenlevel structures similar to that shown in Fig. Sa). 
Here the energy difference of the two paths was controlled 
by fmod, which in effect changes the length of the inter¬ 
ferometer paths in g-space and thus the phase difference 


accumulated by the BEC components which travel along 
i?_|_ and E_. Fig. (b) shows the measured spin po¬ 
larization of the BEC after it has passed both avoided 
crossings, and labeled are the calculated phase difference 
reaching 27r, 47r and Gtt. Similar experiments were run 
for two values of Hoi and the diagram Fig. [^a) shows 
that for smaller Hq, there is a greater energy separations 
of the two eigenlevels as observed in the Stueckelberg 
fringes dependence on fmod- 
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FIG. 2. (a) Two representative modulation induced Stt spin- 
orbit eigenlevels used to measure Stueckelberg interference 
with flo = 1.4Er (blue dashed line) and IJE^ (black solid 
line), and S = 0. The eigenlevels are calculated using 
flc = 0-3Er and fmod = 8 kHz for both, (b) Measured 
Stueckelberg interference fringes for Ho = 1.44?^ (blue) and 
l.lEr (black), with Hm = 0.7Er, and O.SEr respectively 
and ap = IbSOfcr/s for both. For the calculated curves, 
flc = O.SEr, (Ja = 0.07aF, fnp = 0.4, and other parame¬ 
ters match that of experiment. 

In our experiments, a fraction of the ultracold atoms, 
denoted by fnp, do not participate in the Stueckelberg in¬ 
terference due to non-adiabatic initial state preparation. 
In addition, since the BEC is released from its confin¬ 
ing potential at the start of transport it experiences a 
non-uniform acceleration due to atom-atom interactions 
about the central acceleration, ap, created by the ap¬ 
plied force. The non-uniform acceleration is modeled by 
assuming a Gaussian BEC acceleration distribution of 
n{a) = exp [—(a — aF)^/(2(T^)]. The values used 

in this paper for CTq, are consistent with numerically cal¬ 
culated solutions of the Gross-Pitaevski equation using 
a variational method with Gaussian ansatz and param¬ 
eters similar to these experiments |52|. Accounting for 
n(a) and the non-participating fraction, the total spin 
polarization is calculated: 

Stot = (1 - fnp) J n{a)S[a)da (4) 

Including both these effects, we obtain excellent agree- 
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ment with the experiment. This agreement is particu¬ 
larly notable given the period of the modulation is ap¬ 
proximately only an order of magnitude shorter than the 
duration of the experiment, but verifies a time-averaged, 
dressed-state analysis for the eigenlevels can be used to 
model our experiments. 



FIG. 3. (a) Observed Stueckelberg interference fringes for 

different Hm. Data (shown in the order of blue squares, 
red triangles and black circles) were taken with flo = 
1.26,1.40,1.26F;r, S = OEr, Q.m = 0.34,0.70,1.2F;,., and 
ap = 1680fcr/s. Theory fits are shown with a a = O.OTof, 
VLc = 0.14,0.31,0.525,., /„p = 0.4,0.45,0.5. (b) Measured 
fringe contrast M versus modulation amplitude 0.m- Chang¬ 
ing Q.M controls the probability to make a Landau-Zener tran¬ 
sition between E- and E+ , and thus fraction of the BEG that 
splits into each leg of the interferometer. Maximum contrast 
is expected when Plz = 0.5. The theory curve was gener¬ 
ated assuming a linear relation of flc = S2 m/2.25, which was 
fit empirically from the data. For the theory, flo = 1-45^, 
fnp ~ 0.4, <Ja = Q-Qlap- Error bars indicate numerical fitting 
uncertainty of the fringes. 

The gap size, ilc at the pair of avoided crossings af¬ 
fects the probability to travel along either eigenlevel, and 
thus controls the contrast of the Stueckelberg interfer¬ 
ence. Fig. 1^ (a) shows Stueckelberg interference for a 
few representative values of Hm- Defining contrast as 
M = {Smax — Smin)/‘2, Fig. [^(b) shows the contrast ob¬ 
served for various ^Im ■ The results show how ^Im can be 
used to control the Landua-Zener transitions and thus 
the interference fringe amplitude. Given that the max¬ 
imum contrast is obtained when the BEG is coherently 
split with Plz = 0.5, and assuming a linear relation be¬ 
tween ric and Hmj we find ilc = as the best 

fit for these parameters. 

Fig. a shows Stueckelberg interference where the ac¬ 
celeration Of of the BEG, rather than /mod, is varied. 


A secondary optical dipole trapping force was used in 
combination with gravity to control the total accelera¬ 
tion [23] of the BEG and vary the time to transverse the 
two energy paths. Accelerations at which the calculated 
phase accumulations reach 47r and Gtt are labeled. With 
the smaller acceleration, the effect of broadening in ac¬ 
celeration, CTct, due to interactions reduces the observed 
fringe contrast near phase accumulation of 67r compared 
to 47r, but the fringes are still apparent and the model 
agrees well with the experimental results. 

In summary, we have experimentally demonstrated us¬ 
ing SO coupled BECs with modulated Raman coupling 
strength as a method to realize Stt SO coupling and engi¬ 
neer an atomic interferometer. Our theoretical analysis 
matches the results for measured Stueckelberg interfer¬ 
ence fringes, and confirms the treatment of the periodic 
modulation as “dressing” the SO coupled BEG states. 
Interestingly, since the SO coupling is itself the result of 
dressing the mp free-particle dispersion with a Raman 
coupling, this can be considered as “dressing” dressed 
states. These initial experiments show the promise of 
Stt so coupling, which could offer new opportunities 
for studying ground state and phase transitions in a 3- 
minima SO coupled BEG [S] |24| . These also demonstrate 
a different way from other recent studies involving peri¬ 
odic modulation of Raman coupling |25l 126] to engineer 
novel light-induced synthetic gauge fields. 


acceleration (m/s^) 



FIG. 4. Stueckelberg fringes in which the acceleration of is 
varied. The fringe contrast is strongly reduced at smaller ap 
as the relative effect of Ca gets larger with longer total time. 
Do = lAOEr, fmod = 8.5 kHz, Dm = 0.715r, <5 = 0, and 
for theory plots Dc = 0.335^, (Tq = 0.07 x (1680 kr/s), and 
fnp = 0.3. 
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Appendix A: Interferometry theory 

We use the matrix method to solve for the BEG eigen- 
level population resulting from the BEG splitting, phase 
accumulation, and recombination. |'0±) indicates the 
wavefunction in the E± eigenlevels respectively, so the 
state of the BEG is expressed \ijj) = c+\tj}+) -I- c_ |■)/’-)■ 
In operator notation, the state of the BEG is expressed 


The final state after the beam splitter A, phase operator, 
and beam splitter B is thus \tjjf) = Bb^Ba IV’i)- The 
spin polarization is read (V'/l *5 |^), and when solved 


with \tpi) 


0 

1 


results in Eqn. 
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W = (”+), (5) 

and the beam splitters take the form 


—VI — Plz 

\/Plz '\ 

(6) 

VPlz 

Vl— Plz j 

VI — Plz 

VPlz \ 

(7) 

\/Plz 

-VI — Plz ) 


in which P^z is the probability to make a diabatic tran¬ 
sition in the modulation-dressed eigenlevels across the 
avoided crossing, and the negative signs on the diagonals 
account for phase shifts on the wavefunctions at each 
beam splitter [27] . The phase difference accumulated by 
the components of the BEG can be accounted for by a 
phase operator defined: 




e *<^/2 0 \ 

0 j 


( 8 ) 


where (j) is the phase difference accumulated. Readout 
of the final state composition is done by Stern-Gerlach 
separation of the bare-1mi?’) states when the BEG has 
crossed both A and R at a point when the E± eigen¬ 
levels match the bare states to better than 97%, so that 
the spin polarization= - Nmp=-i)/iNmF=o + 

Nm,=-i) « (fV|+) - A^|-))/(A^|+) + iV|-)). Thus, the 
readout of the spin polarization is given by 



FIG. 5. Modulation dressed eigenlevel structure which pos¬ 
sesses a Stt rotation of the spin polarization in the eigenlevels. 
Parameters are the same as used in the experiment shown in 
Fig[^(c), with Uo = 1.3Er, Uc = 0.58i?r, fmod = 10.56 kHz, 
and 5 = 0. 
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